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Abstract. We present the open charm spectra at mid-rapidity from direct 

■ reconstruction of D°, D* and in d+Au collisions at ^/sjvat=200 GeV using the 
, STAR detector at RHIC. The indirect electron/positron measurements via charm 

1^ ' semileptonic decays in p+p and d+Au collisions are also reported. The total cc 

■ cross section per nucleon-nucleon collision is extracted from both direct and indirect 
measurements and are consistent with each other. By combining the D° and 
semileptonic measurements together, the cross section of 1.4±0.2±0.4 mb is higher 

qq . than expectations from PYTHIA and other pQCD calculations. The open charm px 

' distribution from direct measurements covers the pt range up to ~10 GeV/c and 

■ follows a power-law distribution. 
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1. Introduction 

Due to the relatively large mass of the heavy quarks which requires large energies 
for their creaction, the study of heavy flavor productions provides a unique tool to 
test the perturbative Quantum Chromodynamics (pQCD) calculations for heavy quark 
predictions [Tj. The energy dependence of the total cc cross section measured from 
experiments is an important variable to constrain various parameters in the pQCD 
calculations [2 j. It is predicted that the charm quarks are mostly produced from 
the initial fusion of gluons in a high energy hadron-hadron collision jHj. Thus, the 
measurement of open charm hadron production is an important tool to study the initial 
parton distribution function in the colliding nucleons. In relativistic heavy ion collisions, 
the open charm p? spectrum may be modified due to parton energy loss in the deconfined 
hot and dense medium. However, this energy loss is predicted to be significantly smaller 
compared to light quark hadrons due to the "deadcone effect" ptj. As a consequence, an 
increased D/tt ratio as a function of px is expected in central A+A collisions compared 
to p+p collisions at the same beam energy. Therefore, the measurement of the open 
charm px spectrum can provide insights into the properties of the matter formed in the 
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relativistic heavy ion collisions. Furthermore, the production of open charm hadrons in 
the same collision system provides a comparison baseline to test the J /if) production 
mechanism in the relativistic heavy ion collisions at RHIC energies since some theories 
predict a J /if) enhancement due to the charm quark coalescence effect [21 El HI El M while 
others predict a J /if) suppression due to the color screening effect [Til] . 

2. Measurements 

The data used for the direct D°, D* and D ± reconstruction and the charm semileptonic 
decay analysis were taken during the 2003 RHIC run in d+Au and p+p collisions at 
v /J/v r /v=200 GeV in the STAR experiment. A minimum bias d+Au collision trigger was 
defined by requiring at least one spectator neutron in the Au beam outgoing direction 
depositing energy in the Zero Degree Calorimeter (ZDC). A total of 15 million minimum 
bias triggered d+Au collision events were used in the D°, D* and analysis. The 
data samples used in the electron analysis in d+Au and p+p collisions were described 
in Ref. [TT]- 

2.1. Direct D Meson Reconstruction 

The primary tracking device of the STAR detector is the Time Projection Chamber 
(TPC) which measures the kinematics of charged particles. Via the energy loss (dE/dx) 
in the TPC, charged kaons and pions can be identified with their momenta up to 
~0.75 GeV/c. Due to the relatively small decay lengths (cr ~100yum) of the D°, 
D* and D , the limited TPC track projection resolution does not allow the decay 
topology identification by displaced vertices for open charm mesons from the primary 
collision vertex. Thus the event-mixing technique was used to reconstruct the invariant 
mass spectra of the D mesons. Detailed explanation of this method was addressed in 

Ref. [El mad]. 

The low pt {pt <3 GeV/c) D° meson was reconstructed through the decay of 
D° — > K~tt + (D° — > K + tt~) which has a branching ratio of 3.83%. The oppositely 
charged Kit invariant mass distribution after the combinatorial background subtraction 
from event-mixing and a linear residual background subtraction is shown in Panel (a) of 
Fig. with a clear D° signal visible. This distribution was fit with a Gaussian function 
and the mass and width of the D° was found to be 1863+3 MeV/c 2 (1864.5+0.5 MeV/c 2 
in the Particle Data Group [TH]) and 13.8+2.8 MeV/c 2 , respectively. 

The D*^ mesons were reconstructed through the decay of D* + — > D°tt+ (D*~ — > 
D°7T~) with a branching ratio of 68%. Specific treatment for the soft pion daughter 
was performed which required its pt and p between 0.1 to 1 GeV/c and the ratio of the 
D° to 7r s momentum p(D )/p(ir s ) >9 since the n s has an average momentum of about 
50 MeV/c. The invariant mass of the kaon and pion candidates M(Ktt) was calculated 
and required to satisfy the D° mass 1.82< M(Ktt) <1.9 GeV/c 2 . Then a soft pion 
candidate with its charge opposite to that of the kaon candidate was combined with the 
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D° candidate to calculate the invariant mass of a D* candidate M(Ktttt s ). Panel (b) of 
Fig.^shows the distribution of the invariant mass difference AM = M(K7T7T S ) —M(Ktt) 
after background subtractions. A clear signal is seen around the nominal value of 
M(D*) — M(D°). This distribution was fit with a Gaussian function and the mass 
difference and the width was found to be 146.37±0.15 MeV/c 2 (145.42±0.05 MeV/c 2 in 
the PDG) and 0.57±0.16 MeV/c 2 , respectively. 
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Figure 1. The invariant mass distribution after the mixed-event background 
subtraction and a further linear background subtraction of M(Kn) for the low px 
(p T <3 GeV/c) D° + D° in Panel (a), AM = M{K-kti s ) - M(Kn) for the D* ± in 
Panel (b), M(Ktttt) for the D ± in Panel (c), M(Knp°) for the high p T (8.6< p T <12 
GeV/c) D° + D° in Panel (d), respectively. 



The D ± signal was reconstructed through the decay of D + — > i^ _ 7r + 7r + — > 
K + 7t~tt~) with a branching ratio of 9.1%. The invariant mass distribution of M(Ktt7t) 
after background subtractions is shown in Panel (c) of Fig. Q with the D ± rapidity 
\y\ <0.75 and 7.5< px <H GeV/c. A Gaussian function was fit to this distribution and 
the mass and width was found to be 1880±5 MeV/c 2 (1869.3±0.5 MeV/c 2 in the PDG) 
and 16.2±8 MeV/c 2 , respectively. 

The high p T (8.6< px <12 GeV/c) D° was independently reconstructed through 
the decay of D° — * K~7r + p° (D° — > K + 7T~p°) which has a branching ratio 6.2%. The 
analysis was similar to that for the low -D except that an additional 7r + 7r~ pair 
with its invariant mass 0.62< M(tt + tt~) <0.86 was combined with the selected kaon 
and pion candidates for a D° candidate. The invariant mass distribution for M(K7rp°) 
(=M(K-n-K7t)) after background subtractions is shown in Panel (d) of Fig. [TJ This 
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distribution was fit with a Gaussian function and the mass and the width was found to 
be 1850±5 MeV/c 2 and 13.6±4 MeV/c 2 , respectively. 

2.2. Single Electron Analysis 

A prototype time-of-flight system (TOFr) jTJ)J HI] based on the multi-gap resistive plate 
chamber technology was installed in STAR with an azimuthal angle coverage A<p ~ n/30 
and pseudorapidity range -1< rj <0. Besides its capability of hadron identification, 
electrons/positrons could be identified at pt <3 GeV/c by the combination of velocity 
{(3) from TOFr and dE/dx from TPC measurements. In addition, electrons can also 
be identified with 2< p? <4 GeV/c in TPC since hadrons have lower dE/dx due to 
the relativistic rise of the dE/dx of electrons. The inclusive spectra for the low pr 
electrons/positrons measured from TOFr+TPC are shown in Fig. |21 as solid symbols 
and the spectra for higher p? electrons from TPC only are shown as open circles for 
p+p (left) and d+Au (right), respectively. 




p T (GeV/c) p T (GeV/c) 

Figure 2. Upper panels: Electron distributions from p+p (left) and d+Au (right) 
collisions. Bottom panels: Ratios of inclusive electrons over the total backgrounds. 
The gray bands represent the systematic uncertainties. 

The dominant sources of the electron background are the photon conversions 
7 — > e + e~ (dashed curves in Fig.EJ) and 7T° — > r ye + e" Dalitz decay (dotted curves) and rj 
Dalitz decay (dash-dotted curves). To measure these photonic electron background 
spectra, the invariant mass of the e + e~ pairs were constructed from an electron 
(positron) in TOFr and every positron (electron) candidate in TPC. The sum of these 
photonic background is shown as the solid curves in Fig. El In the bottom panel of 
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Fig. |2l the ratio of the inclusive electron/positron spectra and the background is shown 
and clear signal excesses are visible with p? >1 GeV/c. 

3. Results 

3.1. Total cc Cross Section 

From the direct low p? D° reconstruction in the d+Au collisions, the invariant yield 
d 2 N /2npTdpTdy as a function of px after efficiency and acceptance correction was 
extracted in four p? bins at pt <3 GeV/c. Using an exponential fit to the invariant 
yield in transverse mass {raj), the midrapidity yield dN/dy for D° was found to be 
0.028±0.004(stat.)±0.008(s?/s.). We also performed a fit with the combined results of 
D° and electron distributions in d+Au collisions, assuming that the D° spectrum follows 
a power law in pt and that the remaining electrons after the background subtractions are 
charm semileptonic decays. The yield difference between the above two fitting methods 
is much smaller than the statistical uncertainties. 
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Figure 3. The total cc cross section per nucleon-nucleon collision vs. the collision 
energy. The dotted line is a power-law fit. The dashed and dash-dotted curves 
are PYTHIA calculations with different options. The solid curve is a NLO pQCD 
calculation. 

We used the ratio R = N(D°)/N(cc)=0.54±0.05 from e + e~ collisions [15 j to convert 
the D° yield to total cc yield. The d+Au number of binary collisions Nb in and the p+p 
inelastic scattering cross section was used to convert the dN cc /dy in d+Au collisions in 
to da cc /dy in p+p collisions. A factor of 4.7+0.7 jUllin] was used to convert the da/dy 
at midrapidity to the total cross section. The total charm cross section per nucleon- 
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nucleoli interaction for d+Au collisions at 200 GeV is 1.3±0.2±0.4 mb from D° alone 
and 1.4±0.2±0.4 mb from the combined fit of D° and electrons. 

The beam energy dependence of the charm cross section from this analysis is 
depicted in Fig. H by the star symbol and compared to PHENIX [20], UA2 j3T], FNAL 
and SPS j22j measurements. The dotted line is a power law fit, a^ 5 N oc (yfs) n , to the 
data points with n=1.9±0.2, while n ~0.3(0.5) had been observed for charged hadron 
(pion) productions [7] . This indicates a harder behavior of the underlying process going 
from charged hadron production to charm dominated processes. The dashed and dash- 
dotted curves depict PYTHIA calculations with and without higher order processes [T9*] . 
The solid curve depicts the next-to-leading order (NLO) pQCD calculation from Ref. 
[TH] . At a/s=200 GeV, both calculations underpredict the total charm cross section by 
at least a factor of 3. 

3. 2. Open Charm p? Spectrum 

The invariant yield distributions as a function of px for D* ± , and the high p? 
D° were obtained using the same methods as that for the low px D°. With all 
the data points from the direct open charm measurements shown in Fig. 01 the 
Pt distribution was fit to a power-law function, A(l + Pt/Po)' 11 , with the ratio of 
D* / D° = D + / DO as a free parameter, where A, po and n are fit parameters. From the fit, 
we obtained <iiV/c/?/(D )=0.027±0.004±0.007 which is consistent with the exponential 
fit results to the low px D° data points only, (pr)=l-32±0.08±0.16 GeV/c, and the 




Figure 4. The measured invariant yield distributions for D°, _D* ± and and a 
power-law fit to the data points. The D* and data points were scaled by the ratio 
D*/D° = D+/D° = 0.40 ± 0.09 ± 0.13 after the power-law fit. 
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ratio D*/D° = D + / D° = 0.40 ± 0.09 ± 0.13 which agrees with the ratios from previous 
measurements and theoretical predictions [21] within the experimental uncertainties. 



4. Conclusion 

The direct open charm D°, D*^ and signals were reconstructed in d+Au collisions 
at v /J/v r /v=200 GeV using the STAR TPC. Single electron/positron spectra were mea- 
sured using the TOFr+TPC for px <3 GeV/c and using the TPC only for 2< pt <4 
GeV/c in d+Au and p+p collisions at the same beam energy. From the low pt D° and 
the single electron measurements, a total cc cross section was obtained and is at least 
about a factor of 3 larger than PYTHIA and NLO pQCD calculations. The D°, D*± 
and D ± open charm px spectra in d+Au collisions was fit to a power-law fit and the 
D* I D° = D + 1 D° = 0.40 ± 0.09 ± 0.13 was extracted and agrees with the ratios from 
previous measurements and theoretical predictions within the experimental uncertain- 
ties. 
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